Qubit relaxation from evanescent-wave Johnson noise 
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In many quantum computer architectures, the qubits are in close proximity to metallic device 
elements. Metals have a high density of photon modes, and the fields spill out of the bulk metal 
because of the evanescent-wave component. Thus thermal and quantum electromagnetic Johnson- 
type noise from metallic device elements can decohere nearby qubits. In this paper we use quantum 
electrodynamics to compute the strength of this evanescent-wave Johnson noise as a function of 
distance z from a metallic half-space. Previous treatments have shown unphysical divergences at 
z — 0. We remedy this by using a proper non-local dielectric function. Decoherence rates of local 
qubits are proportional to the magnitude of electric or magnetic correlation functions evaluated at 
the qubit position. We present formulas for the decoherence rates. These formulas serve as an 
important constraint on future device architectures. 

PACS numbers: 03.67.-a, 03.65.Yz, 42.50.Lc, 73.21.-b 



Qubits with long relaxation times are necessary for 
quantum computation. Most such devices are controlled 
electrically. This creates a control - isolation dilemma: 
connections from the outside world are what make the 
devices useful, but they are also sources of decoherence. 
In particular, one may wish to place charge or spin qubits 
close to metallic device elements used to confine or con- 
trol the qubits. 

The relaxation of a charge or spin qubit can be induced 
by the thermal and quantum fluctuations of electromag- 
netic fields. The fluctuations of the electromagnetic fields 
are greatly enhanced in the vicinity of conductors be- 
cause of the evanescent waves [lrlal- This evanescent- 
wave Johnson noise (EWJN) has been shown, both the- 
oretically [6[ and experimentally [7[, to be an important 
source of decoherence for atomic qubits near the metallic 
walls of a trap. In this work we investigate the effects 
of metallic device elements in solid-state qubit architec- 
tures. Similar investigations have been carried out pre- 
viously using lumped - circuit calculations of Johnson 
noise |9l4ll|. Here we do the noise calculations taking 
into account the detailed spatial dependence of the fields 
and the important effects of non-local corrections to the 
electromagnetic response functions [8|, [l2[ . 

We pause here to mention that the EWJN originates 
from properties of the metal near its surface. The expres- 
sions for the electromagnetic field fluctuations presented 
in this paper are for a conducting half space. However, we 
have also derived expressions for the strength of EWJN 
in the vicinity of a conducting slab of finite thickness. In 
this case, the electromagnetic fluctuations are indepen- 
dent of the thickness of the slab as long as this thickness 
is significantly larger than the skin depth of the metal. 
As such, we anticipate that EWJN may be alternatively 
interpreted as arising from overdamped surface plasmon 
excititations which exist within a skin depth of the sur- 
face of the metal. |5| 

To describe the decoherence of qubits resulting from 
the evanescent electromagnetic fields surrounding a con- 



ducting gate, it is necessary to compute spectral densities 
of the electric and magnetic field fluctuations. Fermi's 
Golden Rule and the fluctuation-dissipation theorem im- 
ply that the relaxation rate 1/Ti of a (charge or spin, 
respectively) qubit transition of a particular frequency 
will be proportional to the spectral densities at that fre- 
quency. Specifically, the relaxation time T\ of a charge 
qubit with dipole moment d pointing in the ith direction 
at position r and level separation ujz will be given by 
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and Ti of a spin qubit with magnetic dipole moment /2 
in the ith direction at position r and level separation u>z 
will be given by 
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where \u (^ ^ u z ) are the electric and magnetic spec- 
tral densities, respectively. We deal first with fields that 
have been averaged over distances of the order of a and I, 
where a is an interatomic distance and I is the mean free 
path in the electrode. As a result the dielectric function 
e(r , u>) is a local function of space. This approximation 
breaks down in the near vicinity of the conducting sur- 
face, and the influence of a nonlocal dielectric response 
on the correlation function is addressed below. We work 
in a gauge where the scalar potential <f> = 0, so that for 
harmonic fields we have E = iujAJc. 

The spectral densities can be shown [5|, [l3| to be di- 
rectly related to the imaginary part of the equilibrium 
retarded photon Green's function by the relations 
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where i,j are Cartesian indices that run over x, y, z, and 
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The geometry of a particular problem is expressed 
through the function e(r, to). Equations ([TJ and © as- 
sume the charge or spin qubit can be adequately approx- 
imated as a point dipole. The effect of a qubit with 
an extended spatial distribution will be considered in fu- 
ture work. The task of computing Dy in a particular 
geometry is, in general, a complicated problem in elec- 
trodynamics. In this paper we shall limit ourselves to 
the situation where the separation of the qubit from the 
metal surface is much less than any radius of curvature 
of the surface so that the surface can be thought of as 
flat. 

Let z be the distance from the surface. The result for 
the spectral density of the electric field for local electro- 
dynamics has been obtained by Henkel et al |6|: 
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\JuP- jc? — p 2 for p < (J 2 /c 2 and q = 
iy/p 2 — lo 2 /c 2 for p > u 2 /c 2 is the ^-component of the 
wavevector, and p is the transverse component. Our no- 
tation follows that of Ford and Weber 121 . 
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are the Fresnel reflection coefficients. The corresponding 
expressions for the spectral densities of the magnetic field 
are identical to Eqs. ([5]) if we multiply by 1/c 2 and make 
the replacement r 3 O r p . 

We are interested in separations that are sufficiently 
small so that retardation, and hence, radiation of the 
electromagnetic field may be neglected. This is known 
as the quasistatic approximation, and it is formally em- 
ployed by taking the limit c — > oo. This results in the 
greatly simplified expressions 
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FIG. 1: (Color online) Plot of the Ti time of a charge 
qubit computed from Eq. (1) using the local approximation 
(blue) and the full nonlocal theory (red). We used the val- 
ues E F = 7eV, uj = 6tt x 10 8 s~\ v = 6vr x 10 12 s _1 , and 
ujp = 1.6 x 10 16 s _1 , appropriate for a copper surface and a 
device operating in the GHz range. The dipole moment is 
taken as d — |e|as, where |e| is minus the charge on the elec- 
tron and as is the Bohr radius. These results are for zero 
temperature. The dashed horizontal line in the left figure 
represents the strength of the electric field fluctuations inside 
the bulk of a uniform metal. 



Enrploying this approximation has eliminated the func- 
tional difference between the transverse and longitudinal 
components of the field fluctuations. The quasistatic ex- 
pressions differ from the exact values by less than 1% 
when z < 5/10, where 8 is the skin depth of the metal. 
For copper near absolute zero, 5 ~ 3/xm. Equations (O 
diverge as z — > 0, but this divergence is not physical: it 
is an artifact of treating the dielectric function as local 
at distances comparable to the interatomic spacing. In 
any differential equation satisfied by the spatial Fourier 
components of D^ , a local dielectric function will be in- 
dependent of the wavevector while a nonlocal one will 
have a nontrivial wavevector dependence. It is conven- 
tional to represent the spatial Fourier components of the 
nonlocal dielectric function as a tensor quantity in the 
form 
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where we have separated the function into its longitudinal 
ei and transverse et components. 

Nonlocality in the dielectric function changes the re- 
flection coefficients. In the quasistatic approximation r p 
is fH 



i_ eo ^ r dK i 

Ti" Jo k 2 ei(k,uj) 

2 P r, 1 

1 + eo— / an — r 

7T Jo k z ei{k,u}) 



(9a) 




FIG. 2: (Color online) Plot of the Ti time vs. frequency uj of a 
charge qubit computed from Eq. (1) at different temperatures 
(solid red at K and dashed red at 2 K) at fixed z — 10A_f. 
E F = 7eV, v = 6tt x 10 12 s~\ uj p = 1.6 x 10 16 s _1 , and 
d = |e|as, as in Fig. 1. 
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to leading nonvanishing order in the quasistatic ap- 
proximation. Here fc 2 = p 2 + k 2 , 



C-l 



(fc,w) = 1 + -pr-5 : — -—7 r^TTl r" ( 10a ) 

v ; fc 2 v F w + ii>fi((u) + iv) kv F ) 



e t (fc,w) = 1- 



w(oj + if) 



/ t ((w + «/)/A:i; F ) (10b) 



ft{x) 



■ lm> + l)/(x - 1) 



(11a) 



|r 2 -|c(x 2 -l)ln(x + l)/(x-l) (lib) 



v is the electron collision frequency, ui p = (Airne 2 /m) 1 ' 2 
is the plasma frequency, and w F is the Fermi velocity. 
Although Eqs. JSJ) are derived assuming locality, it is a 
convenient fact [12| that these equations are also valid in 
the nonlocal regime, as long as the nonlocal form of the 
Fresnel coefficients Eqs. © are used. 

In Figure 1, we present the zero-temperature results 
for the relaxation time T\ from EWJN for a qubit with 
an electric dipole moment of magnitude |e|as, where as 
is the Bohr radius. Both the local and nonlocal results 
are shown. It is seen that the correct nonlocal dielectric 
function eliminates the unphysical divergence of 1/Ti at 
z = 0. For separations z ~ Xf, the differences are very 
significant, while for z > IQXf, the local and nonlocal re- 
sults nearly coincide. The Fermi wavelength is less than 



a nanometer, so the transition from local to nonlocal be- 
havior occurs well within the quasistatic regime. It is 
interesting to note that for the electric field fluctuations 
there is a crossover region where the nonlocal result be- 
comes slightly larger than the local result in the range 
30A F < z < 3000Af (see Fig. 1), in alignment with the 
results of Volokitin et al [4| who showed an enhancement 
of the nonlocal result above the local result. We see that 
at T = and GHz operations, T\ from spontaneous emis- 
sion is of the order of seconds at separations z ~ 30Af- 
These results are directly applicable to atomic qubits, 
but the rate 1/Ti is proportional to the square of the 
dipole moment, so rates for other charge qubits are eas- 
ily deduced. Figure 2 shows that frequency dependence 
of T\ falls off slowly at higher frequencies, but can be 
very strong at low temperatures T < ftw/kB- 

Figure 3 gives the analogous results for magnetic 
EWJN on a spin qubit with a magnetic dipole moment 
of 1 Bohr magneton. Nonlocal corrections are somewhat 
stronger for this case, and persist to larger distances. In- 
terestingly, the falloff with distance of T\ is slower for 
magnetic EWJN than for electric EWJN. However, mag- 
netic relaxation times are typically somewhat larger than 
electric relaxation times. The crossover of the local and 
nonlocal results is not present in the magnetic case. Fig- 
ure 4 shows that the frequency and temperature depen- 
dence of magnetic EWJN is similar to the electric EWJN 
shown in Fig. 2. Brief mention should be made of the 
dip that is observed as oj — > in Fig. 4. The reflection 
coefficients r s and r p both contribute to the magnetic 
field fluctuations to the same order in lo/c in the qua- 
sistatic approximation. This contrasts with the electric 
case, where only r p contributes to leading order in w/c. 
The dip is a result of competition between the contri- 
butions of r s and r p to the field fluctuations. The r s 
term in xfi is linear in ui as w — > and negative, while 
the r p term is cubic in w as w -> and positive. Wc 
should also mention that for extremely small to the fac- 
tor of coth(Su;/2fcBT) cancels out all u> dependence and 
the dip flattens out as u> — » (not observable in the res- 
olution of Fig. 4) . 

The limit of the nonlocal quasistatic field fluctuations 
as z — > should be of the same order of magnitude as the 
value of these fluctuations inside the metal. To check this, 
we calculate the electromagnetic Green's function inside 
the bulk of a uniform metal using a nonlocal dielectric 
function. The result is 
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FIG. 3: (Color online) Plot of Ti time for a spin qubit at 
zero temperature computed from Eq. (2) in the local approx- 
imation (blue) and the full nonlocal theory (red). Ef = 7eV, 
to = 6tt x 10 8 s~\ v = 6tt x 10 12 s~\ and to p = 1.6 x 10 16 s"\ 
appropriate for a copper surface and a device operating in the 
GHz range. We have taken u = fiB, appropriate for a single 
electron. The rate 1/Ti is proportional to the square of jx, so 
rates for other local magnetic qubits can be easily deduced. 
The dashed horizontal line in the left figure represents the 
strength of the magnetic field fluctuations inside the bulk of 
a uniform metal. 
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We see that there are three relevant distance regimes. 
For z < 30Af, a quasistatic approximation to Henkcl's 
results ((5|) using the nonlocal expression for r p will ac- 
curately describe the field fluctuations. For intermediate 
distances 30\f < z < 5/10, there is a slight enhance- 
ment in the electric field fluctuations from the nonlocal 
expression for r p compared to the local expression. For 
distances z > 5/10, quasistatic local forms (O will accu- 
rately describe the field fluctuations. 

The density of photon states in a metal is very high ow- 
ing to the large polarizability. For blackbody radiation, 
this high density of states does not matter, since total 
internal reflection reduces the outgoing radiation flux to 
its universal Stefan-Boltzmann value. In contrast, the 
evanescent waves are strongly enhanced and the resul- 
tant electromagnetic noise just outside the surface can be 
intense. This is a concern for quantum devices operating 
close to metallic objects. This paper has concentrated 
on the frequency, temperature, and distance dependence 
of the noise, and on the effects of assuming a local di- 
electric function. We conclude that the effect is signif- 
icant for charge qubits with large dipole moments such 
as double quantum dots. The EWJN relaxation may be 
the limiting decoherence effect in designs which involve 
close proximity to bulk metals. For magnetic qubits the 
effects are smaller. We found that nonlocal effects are 
very important at short distances - indeed, local calcu- 
lations can produce spurious divergences. At distances 
large compared to the Fermi wavelength of the metal, 
local approximations work well. 

We have not considered extended qubits for which the 
off-diagonal function Dn(f,r') at r ^ f ' is required. 
We have treated only relaxational decoherence and have 
ignored the possibility of dephasing. These effects arc 
important for many real devices and can be calculated 
using similar methods. 
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FIG. 4: (Color online) Plot of the T\ time vs. frequency of 
spin qubit computed at different temperatures (solid red at 
K and dashed red at 2 K) at fixed z - W\ F . E F = 7eV, 
v = 6ty x 10 1 s~ , and uj p = 1.6 x 10 1 s" 1 , appropriate for a 
copper surface, /i = ub as in Fig. 3. 



Numerical evaluation of (|13p when r = r' gives D xx = 
D zz ~ 3.2 x 10~ 15 Js/m. An evaluation of the Green's 
function outside the metal in the nonlocal quasistatic 
regime for z — > gives D xx ~ 1.32 x 10~ 15 Js/m and 
D zz ~ 2.6 x 10~ 15 Js/m, slightly less than in the bulk of 
a uniform metal, as expected. 
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